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Abstract 

We study transport through one or two ultrasmall quantum dots with discrete 
energy levels to which a time-dependent field is applied (e.g., microwaves). 
The AC field causes photon-assisted tunneling and also transitions between 
discrete energy levels of the dot. We treat the problem by introducing a 
generalization of the rotating-wave approximation to arbitrarily many levels. 
We calculate the dc-current through one dot and find satisfactory agreement 
with recent experiments by Oosterkamp et al. [Phys. Rev. Lett. 78, 1536 
(1997)]. In addition, we propose a novel electron pump consisting of two 
serially coupled single-level quantum dots with a time-dependent interdot 
barrier. 
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I. INTRODUCTION 



Transport through small quantum dots has attracted considerable interest over the last 
couple of years. These quantum dots, small structures formed in a two-dimensional electron 
gas by applying appropriate gate voltages, are characterized by small capacitances to the 
substrate and to the leads connecting them to external voltage sources. Hence there is a 
sizeable charging energy Ec = e 2 /(2C) that has to be provided if electrons are to tunnel 
from the leads to the dot. Transport is blocked at small voltages, a phenomenon dubbed 
the Coulomb blockade since it is a direct consequence of the Coulomb interaction and the 
geometry of the dot. Another aspect that comes up for semiconductor quantum dots as 
opposed to small metallic islands is their discrete single-particle spectrum caused by size 
quantization. 

Many aspects of the Coulomb blockade are now well understood. Recently, a new issue 
has come up, viz., time-dependent transport through small quantum dots. High-frequency 
AC voltages can be applied to mesoscopic structures (e.g., in the form of microwaves). They 
lead to photon-assisted tunneling, i.e., electrons can overcome the Coulomb blockade by 
absorbing photons from the external field. This has become a very active area recently both 
experiment allyiHii and theoreticallyiHl^. 

In this work, we will study transport through an ultrasmall quantum dot with discrete 
energy levels to which a time-dependent field is applied. The electron interaction in the 
dots is taken into account by the Coulomb blockade model. The dots are weakly coupled 
to source and drain reservoirs by tunnel junctions. Time-dependentgate voltages lead 
to photon-assisted tunneling. In contrast to earlier theoretical workBll3 we also take into 
account transitions between discrete energy levels of the dot. In addition, we propose a novel 
electron pump consisting of two serially coupled single-level quantum dots strongly coupled 
by time-dependent fields. 

The paper is organized as follows: in Section |TJ we introduce the Hamiltonian of a sin- 
gle interacting quantum dot with a time-dependent field, connected by tunnel junctions to 
source and drain reservoirs. We discuss the model and its solution by introducing a gen- 
eralized version of the rotating- wave approximation (RWA). In the following section, we 
describe briefly the master equation technique we use to calculate the transport current. In 
this paper, the tunneling is always taken into account by performing a first-order pertur- 
bation expansion in the tunneling matrix element. This is equivalent to consider sequential 
tunneling, assuming the dot to be weakly coupled to the reservoirs such that higher-order 
tunneling processes can be neglected. In Section [TV] we describe the case of a dot with 
two discrete energy levels, which can be solved analytically. Our results for the current are 
presented in Section [V] and compared with recent experiment s3. 

In Section |V| we present a double-dot electron pump, which uses a time-dependent 
interdot barrier as the pumping mechanism. We use the Floquet-matrix technique to find 
a numerical solution to the problem, valid even in situations in which the RWA is not 
applicable. 
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II. MODEL 



As a model for an interacting quantum dot in a time-dependent periodic field cou- 
pled to two reservoirs by tunnel junctions we will use the time-dependent Hamiltonian 
H(t) = H res (t) + H dot (t) + H tu $. Here 

H r es(t) = ]T € ka(T (t) cl aa C ka(T (f) 
k,a,a 

describes noninteracting electrons in the reservoirs {a} = {L,R}, c kaa /c kaa are the cre- 
ation/annihilation operators of an electron with momentum k and spin a in the reservoir a. 
The energies e kaa (t) = e kaa + A a cos(u;i) include a time-dependent shift of the Fermi energy 
of the electrons in the leads due to the applied periodic field, and A a denotes its coupling 
strength to the reservoir a. The Hamiltonian for the interacting electrons in the dot is given 
by 

m<l 

H dot (t) = J2 eiM AA, + E w im{t) Ajma + h.c. + H ch (N dot ) , (2) 

I, it l,m,a 

where AoAio create/annihilate electrons with spin a occupying level / of the discrete, 
equidistant energy levels with level spacing Ae in the dot (I = 1,2, ...,N for a quantum 
dot with N levels). The energy of level I is given by e la (t) = e° CT + A D cos(ut), where the 
time dependence is taken into account by a periodic shift of the level. The coupling strength 
of the field to the dot is given by Ap. The time-dependent transition matrix elements 
w lm (t) = Ai m cos(a>t) describe transitions from level I to level m, i.e., transitions that do not 
change the number of electrons in the dot. The Coulomb interaction between electrons in 
the dot is taken into account by the Coulomb-blockade model 

H ch (N dot ) = E c N 2 dot . (3) 

Here, N dot = J2i,a- AAicr ^ s ^ ne P ar ticle number in the dot, E c = e 2 /2C is the charging energy 
with C = Cl + Cr + C g . Note that we have already taken into account the time- dependent 
part 2N dot n (t) of the original H™ 9 (N dot , t) = E c [N dot +n (t)} 2 by the energies e la (t)E. There 
eno(t) = ClVl^) + CftVnit) + C g V g (t) is related to the polarization charges produced by 
the time-dependent voltages of the left and right reservoirs eVL/R{t) = Hl/r + A L / R cos(u)i) 
as well as the time-dependent gate voltage eV g (t) = n g + A g cos(cjt) applied to the quantum 
dot by the capacitance C g . The tunneling part is given by 

Htun = E T ki cLA* + h.c. , (4) 

k,l,a,<r 

where T kl denotes the tunneling matrix element. 

The time- dependent Schrodinger equation (h = 1) 

i^-m = H \i/>) H (t) = H res (t) + H dot {t) (5) 
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cannot be solved in a closed form due to the time-dependent off-diagonal matrix elements 
w lm (t) in Eq. (g). For 



Ae 



(6) 



one can approximate wi m (t) by 



w lm (t) = A lm cos(ujt) 



A 



W 



I'm 



(t) 



Al m -iuj(l-m)t 



(?) 



This is equivalent to omitting rapidly oscillating terms of frequency uo + Ae, which is much 
larger than uo — Ae as long as (f|) is fulfilled. 

This can be understood as a generalization of the rotating-wave approximation, which 
is well-known in the theory of time-dependent two-level systems (e.g., in nuclear magnetic 
resonance@ or quantum optics^). This generalization can be applied to systems with arbi- 
trarily many levels. It makes it possible to perform a time-dependent unitary transformation 
U(t) = UdV(t) in Eq. (Rl) which removes the time dependence from the non-diagonal ma- 
trix elements of Hdot ty{t)\ an d diagonalizes it afterwards [Ud\. Defining \ip) = U(t)\ip), the 
Schrodinger equation becomes 



with 



ij^) = H (t)\tP) H (t) = H res + H dot (t) 



H^t = UH dot U^ - iU{^) 



(8) 



(9) 



where UH dot W = UH' dot W + H ch (N do t)- The charging part of the Hamiltonian, H ch , stays 
invariant under unitary transformations, since it depends only on the particle number on 
the dot. H' dot is given by the following expression in matrix notation 



HL(t) 
HUt) 



H'dotit) - 

where the submatrices H dot with spin index a = {f, 1} are given by 



(10) 



Hdotit) 



2e lCT (t) 
A 2 e~ i2ujt 



A\e lu}t 
2e 2(T (t) 



A*e i2uJ 
A\e iwt 
2e 3CT (t) 



Here we have assumed Ai m = A\i_ m \ for simplicity. Then V(t) is 



A^^e 1 ^- 1 ^ \ 

i(JV-2)w 
i(JV-3)w 



A* p 
A* p 



:n) 



(12) 



4 



Jujt/2 



V{t) 



( e -^n o 
o 










g gi3ojt/2 







3 i(2iV-3)u;i/2 



(13) 



The eigenenergies (j = 1,2, ...,N) of H dot {t) are then obtained by numerical diagonal- 
ization. The transformed tunneling part of the Hamiltonian is given by 



H tun (t)= £ T^(t)cl a J ja + h.c. 



(14) 



where d^ a jd^ a are the creation/annihilation operators of an electron with spin a that occupies 



the level with energy and 



^•(t)=EW T ^ e ^" 3M/2 - 



(15) 



In addition, for a dot with only two levels the diagonalization can also be performed ana- 
lytically, providing further insight in the underlying physics (see below in Section |TV| ). 



III. THE MASTER-EQUATION APPROACH 

H dot (t) generates the time-evolution operator U (t,t ) = exp(— i Jl o drH (r)), which is 
needed to calculate the tunneling Hamiltonian H^ n (t) = Uo(to,t)H^ un (t)Uo(t,t Q ) in the 
interaction representation. 

The von Neumann-equation ig{t) = [H (t) + H tun (t), g(t)} that describes the time evo- 
lution of the density matrix is also transformed to the interaction representation and q is 
expanded to first order in the tunneling rate. This leads to a master equation for the oc- 
cupation probabilities! P s {t) = (s\g(t)\s) of the occupation number states \s). The states 
\s) represent the occupation numbers of the energy levels of the diagonalized system. In 
the time-averaged dc-case the master equation can be written as a system of coupled linear 
equations 

£[IV-> S - = , (16) 

s' s" 

which can be solved approximately by a suitable truncation. The rate r s /^ s for a transition 
from state \s') to \s) can be expressed as T s >-> s = rf;_^ s + T^_^ s , where 

/J p2tt/uj roo 

^ S = ~ dt dr Re{ ( 8 \H^(t)\^\H^(t - r)\s) } (17) 

71 JO JO 

is the rate associated with tunneling processes from/to reservoir a. H*® n (t) denotes the 
part of the tunneling Hamiltonian ( |T4"D in the interaction picture that corresponds to that 
reservoir. 
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The dc-current through the junction connecting the dot with reservoir a now can be 
expressed in terms of the occupation probabilities and transition rates as well as Nd ot (s) 
(the particle number on the dot while being in state \s)) as 



il = -e £ r2L.[JW) - N dot (s)} P' s 



(18) 



By this means it is possible to numerically calculate the current as a function of transport 
or gate voltage. However, to really understand the resulting I-V-curves, it is helpful to have 
a closer look at the analytically solvable case with only two energy levels in the dot. 

IV. TWO-LEVEL CASE 



The matrix (ITTf) reduces to 



TTcr 
12 dot 



1 / 2e la (t) A*e iuJt 

2 1 Ae"*" 2e 2a (t) 



(19) 



with e la (t) = + Ad cos(ut). The transformation that renders the non-diagonal elements 
time-independent is given by 



V°{t) 



-iuit/2 q 







(20) 



Applying U (t) to the Schrodinger equation leads to 

/2e 1T 



H. 



dot 



\ 



(*)- 


Vu A* 








\ 


A 


2e 2T (t) - 


u 













2e n (t) - 


hw A* 










A 


2e 2l (t) - 


uj) 



(21) 



Calculating Ud and the new energies ej a {t) is now straightforward, they are given by 



-la ^ 2(T ) 



•l) j ^ + A D cos(cut) 



(22) 



where ft = J (Ae — u) 2 + |A| 2 is the Rabi frequency and Ae = e 2a — e\ a . After calculating Ud 

we can write down the creation/annihilation operators dj a /dj cr for an electron that occupies 
energy level j (j = 1, 2) as well as the corresponding time-dependent tunneling matrix 
elements TgAt). Defining the quantity e = Ae — uj ± Vt for to Ae, these matrix elements 



arc 



IAI 



[T^ie^-T^Ae- 



(23) 



T fc ° 2 (t) 



e 2 + IAI 



[T% 2 ee- i " t ' 2 + T% l A*e 



(24) 
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Writing down the tunneling Hamiltonian with these matrix elements and creation/- 
annihilation operators and inserting them into the master equation then permits us to 
calculate the current in a straightforward way. It turns out that there are two possible 
ways for an electron to tunnel through the dot, corresponding to the two terms with e ±wt l 2 
in the tunneling matrix elements (Fig. [[]). The transport peaks present in the absence of 
a time-dependent field are split in a two-peak group (with peaks at distance u) that also 
shifts its position due to the A-dependence of the energies 6j a (t). 



V. RESULTS 

For a dot with two spin-degenerate levels four groups of current peaks will appear in the 
I-V-curve, separated from each other by 2Ec- If the field- induced inner transitions between 
the levels are neglected (all A\i_ m i = 0), there is one main peak per group accompanied by 
Bessel-type sidebands at separations nu (with n = ±1, ±2, ±3, ...). These side peaks are due 
to photon-assisted tunneling (PAT)i. The existence of these side peaks has been recently 
verified experiment ally0. 

If inner transitions are taken into account, our calculations show that the single main 
peak will shift with increasing An m i. In addition, N — 1 peaks will appear at distances 
nuo (with n = 1, N — 1), see also Fig. [1|. This leads to a totally different picture of the 
current-peak positions and heights, the weight of the peaks shifts as well as their positions. 

To visualize only the influence of inner transitions, we set Ad = 0, i.e, there is no PAT, 
and plot the peak group corresponding to a dot occupied by one additional electron as a 
function of the increasing strength of the inner transitions (increasing Aij_ m i). Two cases, 
one with two and one with three degenerate energy levels are shown (Fig. |3|). The shift of 
the main peak and the appearance of the one (two) additional peak(s) as described above 
is clearly visible. 

Of course it is also possible to include the inner transitions in the master equation 
in a perturbative way. To do this, transition rates between the levels (similar to those 
describing the tunneling of electrons to and from the dot) have to be calculated in first- 
order perturbation theory. Then the peaks do not shift with increasing Au_ m i, because the 
master equation is written in the basis of the unperturbed states. But also the effect on 
the peak heights is significantly smaller than in our modified rotating-wave approximation 
approach. In Fig. |2] this is illustrated by plotting again the peak group for a dot with two 
levels and one extra electron. Again we set Ad = 0, i.e., neglecting PAT. For a fixed, 
relatively small value of A, the peak shift is not significant. 

In recent transport measurements on ultrasmall quantum dots with applied microwaves 
such a peak group was studiecfi A dot with two levels contributing to transport was used. 
The first side peak rises more strongly with increasing microwave power than expected from 
the PAT model. In Fig. [| we compare the measured I-V-curves with our results. For small 
values of A our model is in good agreement with the experiment for both microwave fre- 
quencies. In particular it describes correctly the strong increase of the right side peak. 
The numerical calculations have been done using non-degenerate levels because in the ex- 
periments a strong magnetic field (B = 0.91T) (which was supposed to suppress plasmon 
excitations) lifted the spin degeneracy. 
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VI. DOUBLE QUANTUM-DOT ELECTRON PUMP: FLOQUET-MATRIX 

APPROACH 



The single quantum dot with two discrete energy levels connected by a matrix element 
wimit) = A + A cos (cut) can be mapped to a double dot system where two dots with one 
energy level are strongly coupled with each other. But in the formalism discussed above it 
would not be possible in this case to choose the gate voltages and the microwave frequency 
arbitrarily, due to the restrictions imposed by the rotating-wave approximation. Also, the 
time-dependence of the two gate voltages would have to be the same for both dots. 

In this section we discuss a double dot system with time-dependent gate voltages that 
differ by a relative phase if. The dots are strongly coupled by a time-dependent tunneling 
barrier . Such a system connected to reservoirs as shown in Fig. [| can work as a electron 
pump, resulting in a current even if no transport voltage is applied. We generalize the 
work of Stafford and Wingreen0 to the case where the height of the tunneling barrier also 
depends on time. The time-dependent Schrodinger equation is solved by the Floquet-matrix 
approachB. In general, according to Floquet's theorem, a differential equation with periodic 
coefficients like equation (|5|) has solutions of the form (for a dot with N levels)H 

|Vi(*)>=e-^*|Vi(*)> (J = 1,-,A0, (25) 

which inserted in the Schrodinger equation (|5[) result in an eigenvalue problem for the states 
\<Pi) 

(tf(t)-^)|^) = e>,>. (26) 

The states \ipj) have the same periodicity as the dot Hamiltonian, i.e., \<fj(t + 2it/uj)) = 
\(Pj(t)). Due to this it is possible to expand H dot {t) and \fj(t)) in a Fourier series 

oo 

H*M = E <U iWt (27) 



oo 

\<P&)) = E e im3t \<pf) • (28) 

n=— oo 

In the basis { | / ) } of the eigenstates of the diagonal dot Hamiltonian H® ot = Yji, a e ^Aa^iu 
with uncoupled energy levels e° (i.e., H% ot \ I) = e°| /)), the solution \ipj(t)) can be expressed 
as 

N N oo 

lV>i(*)> = E( * l*>i(*)> ^ 1 1) = E E e-^V-* tpff 1 1 ) . (29) 

1=1 1=1 n=—oo 

If we insert ( |2"T| ) and in the time-dependent Schrodinger equation (|), multiply with the 
bra (i | from the left and define the matrix elements = (i \H^ t \ /), we get an infinite 
system of coupled linear equations describing the eigenvalue problem for the quasi-energies 
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N oo 

E E [Ht k) + ruvSM $ = . (30) 

1=1 k=—oo 



If H dot (t) is given by (0) with wi m (t) = A 0lm + A im cos(ut), this becomes 



N oo 



A; 



E E [(( e ° + nUJ )$il + A 0u)$nk + "y (<5n,fc+l + <W-l)] (f$ = Zj^f . (31) 



1=1 k=— oo 



The quasi-energies e,- and the eigenvector components (p\j can be calculated numerically by 
truncating this infinite system of coupled equations at a sufficiently large finite n. It is now 
possible to transform the tunneling Hamiltonian and calculate the current by the master 
equation technique analogous to the previous sections. The resulting dependence of the 
pumped current on the frequency and the amplitude of the applied microwaves is plotted in 
Fig. H and Fig. [7]. Here and in the rest of the paper we used E ch = 75T, T L = T R = T. In 
Fig. | and Fig. |we set e x = -10T, e 2 = 10r, fi L = fi R = 0, T = 5T. 

Figure |^ illustrates the situation where the microwaves couple only to the interdot barrier. 
There are current peaks if the photon energy equals an integer fraction of the quasi-energy 
level splitting. Because the quasi-energies themselves depend on the amplitude A (analogous 
to the RWA calculations above), the current peaks shift to higher photon energies with 
increasing A. This plot illustrates the effect created purely by a time- dependent barrier. 

However, in a real experiment the microwaves would also couple to the gate electrodes 
and the interdot coupling would have a finite time-independent part, i.e., wi m (t) = A + 
Acos(o;t). This is shown in Fig. 0. It is clearly seen that the overall peak height increases 
compared to Fig. ||. 

The behavior of the current through an interacting double quantum dot for finite trans- 
port voltages and two independently varied gate voltages is plotted in Fig. || and Fig. || 
(Hl = 2.5r, ji R = — 2.5r, T = 3r in both figures). Figure |S] illustrates the case where the 
microwaves are coupled to the interdot barrier only, i.e., with a time-dependent coupling 
matrix element. Instead, in Fig. |S] the microwaves are coupled to the gate electrodes, as- 
suming a static interdot matrix element. As can be clearly seen, the same value for the 
coupling matrix element in both case leads to a significant increase of I max - 

In conclusion, we have calculated the photon-assisted transport current through a single 
interacting quantum dot with an arbitrary number of discrete energy levels. We have taken 
into account field-induced inner transitions in a non-perturbative way by generalizing the 
rotating-wave approximation to more then two energy levels. We compare our results to 
recent experimentsB and provide an explanation for the unexpected height of the first side- 
band current peak. 

We would like to thank J. Konig and G. Schon for discussions and suggestions and 
especially T. H. Oosterkamp and L. P. Kouwenhoven for providing us with the experimental 
data. The support by the Deutsche Forschungsgemeinschaft, through SFB 195, is gratefully 
acknowledged. 
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FIGURES 



7fr 

e 2 i +C0/2 




left reservoir quantum dot right reservoir 

FIG. 1. Energy landscape of a quantum dot with two non-degenerate levels and inner transitions 
induced by a time-dependent field of frequency uj. The effective energy levels ?j are shown (solid 
lines). Note that electron transport occurs, when the chemical potentials Hl/r match one of the 
quasi- levels (dashed lines) shifted from the energies €j by ±uj/2 due to the time-dependent tunneling 
matrix elements T£-(t). 



0.04 



0.03 




FIG. 2. /d c -V^-curve for a dot with two degenerate levels. The RWA model (solid curve) is 
compared with first-order perturbation theory for calculating the transition rates between the two 
energy levels (dashed curve). In RWA the first side peak is significantly enhanced. t\ = — 25F, 

e 2 = 25T, uj = 6or, t = 5r, r 2 = r 2 = o.sr, a d = o, a = lor. 
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"We 

FIG. 3. Evolution of the current-peak group corresponding to one electron occupying a dot with 
two (top) and three (bottom) spin-degenerate levels with increasing A. PAT has been omitted. 
A D = 0, Ae = 95r, uj = 55I\ T = 3I\ 
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FIG. 4. Measured Idc-Vg-curves (solid lines) (Oosterkamp et aZ.S) with applied microwaves 
of frequencies 42 GHz (left column) and 61.45 GHz (right column) versus our theoretical results 
(dashed lines). Left column: Ti = 0.2I\ A D = 1.875A. Right column: Ti = O.ir, A D = 3.5A. 
In both cases, Al = — j^Ad, Ar = -^Afl. Parameters determined by the experiment: 
ei = -47.5r, e 2 = 47.5r, w = 135r (left), oj = 200r (right), and V = 1.287/xeV. 
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FIG. 5. Energy landscape of two serially coupled quantum dots connected by a weak 
time-dependent barrier (periodicity lo). A additional time dependence may also be applied to 
the gate electrodes of both dots. With (il = fJ-R and applied microwaves the system works as an 
electron pump, pumping electrons "uphill" from left to right. 
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FIG. 6. Current response of a double quantum-dot electron pump versus frequency co for a 
time-dependent barrier separating the dots. Aq = 0, Ad = 0. 
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FIG. 7. Same as Fig. ^ with time-dependent gate voltages and assuming a phase difference of 
7r between the two gates. Aq = 2T. 
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FIG. 8. Dependence of the transport current on the positions of the energy levels, i.e., the gate 
voltages, for a time-dependent interdot barrier. Aq = 2r, A = 3r, A^ = 0. 
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FIG. 9. Same as Fig. || for static interdot coupling (A = 0) and time-dependent gate voltages, 
assuming a phase difference of it between the two gates. Aq = 2T, A = 0, Ad = 3T. 
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